Phase transitions in dipolar gases in optical lattices 
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We investigate the phase diagrams of two-dimensional lattice dipole systems with variable geom- 
etry. For bipartite square and triangular lattices with tunable vertical sublattice separation, we find 
rich phase diagrams featuring a sequence of easy-plane magnetically ordered phases separated by 
incommensurate spin-wave states. 
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A recent breakthrough in the cooling of dipolar gases 
in optical lattices [lj , following a decade of intensive re- 
search dHl], opens a door into the earlier inaccessible 
many-body physics of lattice systems with anisotropic 
long-range interaction. Although bulk crystalline dipo- 
lar systems are abundant in nature fiol their experi- 
mental investigation has been hindered by the extremely 
low temperatures required for the observation of ordering 
transitions [l2| and the absence of continuously variable 
parameters. Recently, artificial two-dimensional dipolar 
systems, such as lithographically created nanomagnet ar- 
rays, have been realized (l3lll4j ). However, it is the advent 
of optical lattices with tunable lattice structure contain- 
ing ultracold dipolar gases, that creates numerous possi- 
bilities for studies of previously unexplored phase tran- 
sitions - both classical and quantum - between ordered 
and disordered phases of this fundamental many-body 
system. 

In this Article we analyze a series of magnetic phase 
transitions in a classical dipolar gas in deep optical lat- 
tices [square, Fig. [T] and triangular, Fig. f2 obtained 
from bipartite monolayer lattices by vertical separation, 
z, of the two sublattices. One way to realize such systems 
would be loading of ferromagnetic spinor Bose-Einstein 
mini-condensates in the nodes [i~5| of a deep bilayer op- 
tical lattice created with the help of a painted poten- 
tial technique [ijl], which would allow for high degree of 
control over the shapes of optical lattices and interlayer 
separation. 

We find that, upon the variation of z, each system ex- 
periences a sequence of easy-plane magnetically ordered 
phases separated by incommensurate spin- wave states, 
which could be detected with the help of Bragg diffraction 
of light (l"7l - [l9t . The phase diagram for the square lattice 
on the z — T plane is shown in Fig. [TJ For sufficiently 
small separations z <C a, where a is lattice constant, we 
reproduce the earlier predicted (2fl [2l| canted antiferro- 
magnetic phase, AFk, with the ordering vector K. For 
z > a, we find an antiferromagnetic phase, AFm, with 
larger unit cell and ordering wave vector at the M-point 
of the Brillouin zone of the bipartite lattice. For interme- 
diate interlayer distances, we find a stable ferromagnetic 
phase (F), separated from the antiferromagnetic ones by 
incommensurate spin- wave states (ISW). At the critical 
temperature T c , all of the ordered phases feature a degen- 
eracy in the orientation of magnetization, characterized 
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Figure 1: (color online) (a) Dipolar magnetic gas on a bipar- 
tite square optical lattice as seen along [0,1] axis and the first 
Brillouin zone. The two sublattices, A and B, are vertically 
separated by the distance z. (b) Phase diagram of the dipolar 
gas: three commensurate phases (AFk, F, and AFm) sepa- 
rated by two incommensurate spin-wave phases (ISW1 and 
ISW2) with phase boundaries at z\ = 0.245a, Z2 = 0.262a, 
zz — 0.629a and Z4 = 1.061a, where a is lattice constant, and 
To = /io/i 2 /47ra 3 . 

T c , such a degeneracy is lifted, and Fig. Q] shows the 
optimal orientation of the order parameter for the low- 
temperature states. The structure of the intermediate 
incommensurate phases, ISW1 and ISW2 (Fig. [!} has 
also been established for T —> T c , while their nature at 
low temperatures remains an open question. 

We find that the phase diagram for the bipartite tri- 
angular lattice (which forms a honeycomb lattice when 
z = 0) also contains a sequence of commensurate and 
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Figure 2: (a) Dipolar magnetic gas on a bipartite triangular 
optical lattice as seen along [0,1] axis and the first Brillouin 
zone. The two sublattices, A and B, are vertically shifted 
separated by the distance z (at z = 0, they form a honey- 
comb lattice), (b) Phase diagram of the dipolar gas: three 
commensurate phases (Hk, AFm, and F) separated by two in- 
commesurate spin-wave phases (ISW3 and ISW4) with phase 
boundaries at Z\ — 0.256a, Z2 — 0.2600a, z$ = 0.479a and 
z C 4 = 0.577a, where a is lattice constant and To = fiofJ 2 /47ra 3 . 



incommensurate magnetic phases, FigJ5J For z < a, we 
find a helical phase with the ordering vector K (Hk), 
which is specific for dipoles on a 2D honeycomb lattice 
[20I ] . A large vertical displacement of the two sublattices 
of the honeycomb lattice results in two weakly coupled 
triangular lattices, for which the ground state is ferro- 
magnetic (F) (2^1 - In between those two extremes lies 
an antiferromagnetic phase AFm with the ordering vec- 
tor M [22j , separated from the helical and ferromagnetic 
phases by parametric intervals, where the magnetization 
texture is incommensurate with the lattice. 

To find phase diagrams in Figs. [TJ [21 we consider an 
ensemble of classical magnetic dipoles /is, (|si| 2 = 1), 
placed on the A or B sites r™ (n = A, B) of the square 
(Fig. [JJ or triangular (Fig. [2]) lattices, with interaction 



energy 

H =\H E E J^(r?-rf)s a (r^(vf), 

Here is the vacuum permeability. The Hamiltonian 
is invariant under the group Q of simultaneous rotation 
of magnetic moments in the xy-plane and the lattice ro- 
tations through the angles 7r/4 for square and 7r/3 for 
triangular lattices. 

In order to identify the thermodynamic average, 
(s(r™)) (n = A, B), of magnetization for various inter- 
layer separations z, we apply the Landau theory and 
study the free energy in the vicinity of the transition 
temperature (23| . 

q 

expressed in terms of a 6-vector 

S(q) T = (5I(q), S\(q), S%(a), S%(q), 5|(q), 5|(q)), 

where S" (q) is the Fourier transform of the order param- 
eter (s(r™)) of magnetic moments on A and B sublattices. 
In Eq. [2J 1 is a 6 x 6 unit matrix, and F^ m > incorporates 
higher-order invariants under the group Q built using the 
order parameter (m = 4 for square and m = 6 for trian- 
gular lattices). A 6 x 6 matrix J has elements 

■^(q) = (W E jaP W ~ r?V q - (rr -A (3) 

where N is the number of unit cells. For each wave vec- 
tor q, matrix J has 6 eigenvalues, A, and 6 eigenvec- 
tors, V(q). The lowest negative eigenvalue Ao(qo) found 
among A 7 (q) by varying wave vector q over the Brillouin 
zone determines the polarizations V(qo) and the wave 
vector qo of the most favorable magnetic ordering and 
the transition temperature 

T c = -imin(A 7 (q)) = i|A (qo)|. (4) 

In Fig. [31 we show plots for A 7 (q) for the square bi- 
partite lattice with various vertical A-B sublattice sepa- 
rations. For < z < Zi, where z\ = 0.245a (a is lattice 
constant) qo coincides with one of K-points of the Bril- 
louin zone as shown in Fig[3] (a)-(b). This corresponds 
to the AFk phase (Fig. [TJ) with the order parameter 

(s(r™)) oc (cos(K • r™ + 9),psm(K • r? + 9)), (5) 

where p = ±1 for n = A/B, and 9 is a constant. In Eq. 
. the degeneracy in 9 is lifted by the higher-order terms 
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appearing after taking into account thermal fluctu- 
ations. As z increases from z\ to zi = 0.262a, qo con- 
tinuously moves from K- to T- point (Fig0](a)), and the 
corresponding eigenvector Vo(qo) determines the mag- 
netization texture 



(s(if )) oc z x q cos(q • if) 



(6) 



of the incommensurate phase ISW1 illustrated in Fig. 
[T] (b), where z is the unit vector perpendicular to the 
plane of the lattice and qo = (<Zo,0). For z^ < z < 
Z3, where z% = 0.629a, qo lies at T-point (Fig. [3] (d)), 
which corresponds to the easy-plane ferromagnetic (F) 
ordering. As z increases from z^ to z^ — 1.061a, qo 
continuously moves from T- to M -point (Fig. 0] (b)), 
which determines the incommensurate spin-wave state 
ISW2 and the order parameter given by Eq. © with 
qo = (1/v / 2)(9oj <Zo)- For z > z&, qo is at one of M-points 
(Fig. El (f)) giving rise to the phase AFm, which can be 
viewed as two weakly coupled "AFk" phases on each of 
two square sublattice. The form of the order parameter 
in each of the commensurate phases is given in Table |TJ 
The phase diagram for the bipartite triangular lattice 
(Fig. [U (b) with order parameters listed in Table fl} is 
somewhat similar to that for the square bipartite lattice. 
For < z < zi, where z\ = 0.256a, we find that qo is at 
one of the K-points. This corresponds to the helical Hk 
phase with the order parameter given by Eq. ([5]), where 
vector K is at the corner of the hexagonal Brillouin zone 
of triangular lattice (Fig. [U (a)). Such a phase has been 
discussed in relation with a dipolar gas in a planar honey- 
comb lattice [2(j • For z\ < z < z%, where z-i = 0.260a, qo 
continuously shifts from K to M-point giving rise to the 
incommensurate phase ISW3 with magnetization texture 

( s ( r ™)> « (-sin(qo • r?),-pa v cos(q • r?), a z sin(q • rf )), 

where < a y and a z <C a y are z-dependent, qo = 
g (K(l - c) + Mc)/|(K(1 - c) + Mc)| (0 < c < 1). For 
Z2 < z < Zs, where 23 = 0.479, qo lies at one of the 
M-points, which corresponds to the easy-plane antifer- 
romagnetic phase AFm shown in Fig. [2] (b) [22j |. For 
Z3 < z < Z4, where 2:4 = 0.577a, qo moves from M- to T- 
points, which determines the incommensurate spin-wave 
state ISW4, 



(s(r™)> oc (sin(q • r"),cos(q • r?)) 



(7) 



with qo = qoM/|M|. Finally, for z > Z4, qo lies at T- 
point, which corresponds to an easy-plane ferromagnetic 
state. In the limit z — > 00, this coincides with the ground 
state calculated for a dipolar magnet on a plane triangu- 
lar lattice f20l |. 

The above analysis of magnetic phases of dipolar gases 
on square and triangular bipartite lattices, limited to the 
quadratic terms in the Landau theory, is formally valid 
at T — >• T c . To extend the phase diagrams in Figs. [T]and 
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Figure 3: Eigenvalues, Ai(q), of the dipolar tensor, J™„/, for 
the square lattice as functions of the wave vector q along 
symmetric directions in the Brillouin zone (see Figflja)) for a 
set of lattice displacements, z: (a) 2 = 0; (b) z = 0.140a; (c) 
z = 0.248a; (d) z = 0.500a; (e) z = 0.707a; (f) z = 2.120a. 
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Figure 4: The minimal eigenvalues Ao vs. q along symmetric 
directions in the Brillouin zone for incommensurate phases: 
(a) ISW1 and (b) ISW2 for a representative set of z: (A) z = 
0.245a; (B) z = 0.250a; (C) z = 0.255a; (D) z = 0.260a; (E) 
z = 0.262a; (F) z = 0.629a; (G) z = 0.651a; (H) z = 0.743a; 
(I) z = 0.849a; (J) z = 1.061a. Arrows show the positions of 
the minima, qo, of A. The ordering vector, qo continuously 
moves along the straight lines connecting the K- and F- points 
(a) and T- and M-points (b), as z increases. 



[5] to low temperatures, we investigate the stability of the 
ordering patterns described in Table U near T = using 
the linear spin- wave theory [23j j . That is, we expand the 
interaction energy in small deviations of on-site magnetic 
moments from the ground-state value (s(rf)), 



£ = £(S) + ^5>t(q)A(qMq) 



(8) 



Such deviations have to respect the constraint |s(rf)| 2 
1 and can be parametrized as 

s(r*) = (s(r*)) + (z * 

-(l/2)(s(r?))[((7||(r*)) 2 + ((7,(r?)) 2 ]. 
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Table I: Order parameter for each of the commensurate phases of a dipolar gas on bipartite square and triangular lattices. 



Lattice 


Phase 


S = {s(r?))/|(s(r?))| 


qo 


So 


Square 
n = A,B 


AF K 


(cos(K ■ r™ + 6>),psin(K • r? + 9)), p = ±1 for A/B 


K = (0,^) 


2^, fc € Z 


F 


(cos 9, sin 9) 







AF M 


(cos(M ■ r? ± n + tt/4), sin(M ■ r? ± 6>„ 4- tt/4)) [24J 


M = (^,^) 




Triangular 
n = A,B 


H k 


(cos(K ■ r™ + 0),psin(K • if + 9)) 


K=(0,^) 




AF m 


(sin(M-r"),cos(M ■!■")) 


M = (^,0) 




F 


(cos 9, sin 9) 





2^, fc € Z 



Here we use the Fourier transform, c(q) T = 
(<rj|(q),^(q),...,CT[j(q),<7*(q),...) 01 ^a^i); index fc = 
1, ...,M labels sites within the magnetic unit cell of 
a commensurate phase, so that the Fourier transform 
of (s(r*)) is Sfc(q) = Sfc<5(q) and Nm is the num- 
ber of magnetic unit cells. The 2M x 2M matrix 
A(q) has elements A^'J(q) = £ Q(3 (z x S fc ) Q J^(q)(z x 
S fc ,)" + /(0), A^(q) = £ Q (z x S fe )«J- (q), 4>IJ(q) = 
J^(q)(z x S k ,f, A££,(q) = J&(q) 4- /(0), where 
/(0) = (-1/2) J2 a p S fc a J fe t(0)5f,. We find that all eigen- 
values of A(q) are positive for q 7^ within the same 
intervals Zi < z < Zi+\ , whereas at the edges of the 
intervals, the lowest eigenvalue of A(q) acquires nega- 
tive values reflecting a tendency towards incommensurate 
textures. 

For most of the phases in Table HI the interaction is 
degenerate in angle 9 (Figs. [T]and[5]). This degeneracy is 
lifted [1^1 by thermal fluctuations leading to the higher- 
order expansion terms in the Landau theory. To find, at 
least for T <C T c , such symmetry breaking contributions, 
we take into account fluctuations er(q) of the order pa- 
rameter following [25| and calculate the entropy part of 
the free energy 

8F S = -TlnJ] / IMq)exp[-iat(q)A(q)a(q)] 
q J 

E MdetA(q)). (9) 

We evaluate the dominant contribution from the fluctu- 
ations lifting the degeneracy with respect to angle 9: 

5F S w TA - TA\ cos(m<2), A , A x > 0, (10) 

where m = 4 for square and m = 6 for triangular lattice 
[2^ . This determines the optimal choice 9q shown in 
Table |TJ For T T c , such entropy terms give rise to the 
crystalline anisotropy contribution cx T c |S|" l cosm0 
in Eq. ©. 

For magnetic dipolar gases in deep bipartite (bilayer) 
square and triangular optical lattices, the predicted phase 



diagram may appear very much within the experimen- 
tally accessible range of controlled parameters. For deep 
optical lattices with a ~ 1/im and optical field trapping 
mini-condensates of 10 3 spin-aligned 87 Rb atoms per unit 
cell, we estimate T c ~ 50^K in the phase diagram in Figs. 
Q] and [5J Moreover, as the electric and magnetic dipole 
interactions are mathematically equivalent, the phase di- 
agram in Figs. [T] and [5] should be applicable to the elec- 
tric dipolar systems, where we estimate T c ~ lOOnK for 
ferro- and antiferroelectric transitions in molecules with 
a dipole moment d ^ ID. 
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